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Abstract Clusters of galaxies are the most massive objects in the Universe and precise 
knowledge of their mass structure is important to understand the history of structure 
formation and constrain still unknown types of dark contents of the Universe. X-ray 
spectroscopy of galaxy clusters provides rich information about the physical state of 
hot intracluster gas and the underlying potential structure. In this paper, starting from 
the basic description of clusters under equilibrium conditions, we review properties of 
clusters revealed primarily through X-ray observations considering their thermal and 
dynamical evolutions. The future prospects of cluster studies using upcoming X-ray 
missions are also mentioned. 
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1 INTRODUCTION 

According to the standard cosmological model, the Universe began 13.8 billion years ago, and con- 
sists of 4% bar yonic matter, 23% dark matter (of unkno wn type) and 73% dark energy (also of un- 
known origin) ( Larson et al.ll2.01 ltlKomatsu et al.ll2.01 lb . Through interactions of these constituents, 



the associated cosmic structures have been evolving up to now. Our description of the Universe 
is often based on the notion that large objects, like galaxy clusters, that formed out of the evolv- 
ing large-scale structure, have attained an equilibrium state in their matter and energy constituents. 
However, is this truly a natural assumption? To tackle this problem, by focusing on objects appearing 
at the top of the hierarchical structure formation, namely clusters of galaxies, is vital in astrophysics. 

Clusters of galaxies are the largest gravitationally bound systems in the Universe. This makes 
them very important probes of cosmology. Thus a precise knowledge of their mass structure is very 
important to measure the large-scale structure and to test cosmological models. In visible light, they 
are identified as groups of ~ 100— 1000 galaxies, extending over ~ 10 7 light years (Figure[T]left). On 
the other hand, X-ray observations of clusters have drastically changed our view of cosmic structure: 
hot gas fills inter-galactic space and emits strong X-rays (Figure[TJright). Furthermore, the total mass 
of hot gas exceeds the sum of galaxy mass by two-three times. To confine the hot gas by gravitational 
forces, invisible matter, "dark matter," of five times larger mass is required. As techniques in X-ray 
spectroscopy and imaging observations progressed, the presence of a complex temperature structure 
was also found in the X-ray emitting gas. Those facts have revealed that clusters preserve the past 
history of being built through complex interactions, particularly merging, between smaller systems. 
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Fig. 1 Optical (left) and X-ray (right) images of a cluster of galaxies, RX J1347.5-1 145, 
taken with the Hubble Space Telescope (the Multi-mission Archive at STScI) and the 
Chandra satellite, respectively. In both panels, a side of the figure is 1 10", corresponding 
to about 630 kpc. The white circle in the right panel indicates a location where extremely 
hot thermal gas has been discovered (see Section l4~3l . 



Thus the clusters are no longer thought to be in an equilibrium state, but rather dynamically evolving 
on cosmological time scales. 

This paper is organized as follows: in Section 2 a general description of clusters is summa- 
rized. In Section 3-4, properties of clusters of galaxies revealed primarily by X-ray observations are 
reviewed in light of their thermal and dynamical evolutions. Finally in Section5, future prospects 
are briefly mentioned. We use Qjvf = 0.3, = 0.7 and hjo = Hq/(70 kms -1 Mpc -1 ) = 1 
throughout the paper except where noted. 



2 SPATIAL DISTRIBUTION OF MATTER AND X-RAY EMISSION 

In the X-ray energy band, clusters of galaxies look very different from the optical view; hot diffuse 
plasma with a temperature of ~ 10 7 — 10 8 K fills the intracluster space. The X-ray emitting hot 
plasma is confined in the cluster gravitational potential and is believed to trace the underlying dark 
matter distribution. The typical X-ray luminosity of clusters is 10 44 ~ 45 ergs -1 , and the electron 
number density of hot plasma at the center of clusters is typically 10~ 3 ~~ 2 cm~ 3 . In what follows, 
the general view of clusters of galaxies under equilibrium models is summarized. 



2.1 Hydrostatic equilibrium condition and the /3-model 

Since the collision time scales for ions and electrons in the intra cluster gas are much shorter than the 
time scales of heating or cooling, we can treat the gas as a fluid dSarazi n 1988). In general, the sound 
crossing time t s , i.e., the time required for a sound wave in the intracluster gas to cross a cluster with 
radius R, 

2R „ „ ( R 



ts c s ~ 2Gyl ^iMpcJ V 1000 km s- 1 J (1) 

is shorter than the probable age of the cluster or the Hubble time, tjj = Hq ~ 14 Gyr. Thus the gas 
is considered to be in hydrostatic equilibrium. In addition, if the cluster is spherically distributed, the 
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hydrostatic equation reads 



Pg dr 



d(j> _ GM (r) 
dr r 2 



(2) 



where M(r) is the total cluster mass (i.e., dark matter + galaxies + hot gas) within the radius r 
and P g is the thermal pressure and a product of gas density and temperature, n g {r)kT g (r). If the 
self-gravity of the gas is ignored, the distribution of gas is determined by the cluster potential, <j>. 

A temperature gradient in the plasma is smoothed by heat conduction. If the heat conduction 
were sufficiently rapid compared to other important time scales, the gas would become isothermal. 
Substituting the gas pressure P(r) = n g (r)T g into Eq. (1) and assuming T g is constant, we obtain 



<iln p g 

dr 



Hrrip d(j) 



(3) 



where /i is the mean molecular weight, ~ 0.6. Similarly, the galaxies are bounded in the gravitational 
potential, whose hydrostatic condition is written as 



d\n p G 

dr 



1 d<P 
a 2 dr 



(4) 



a is the line-of-sight velocity dispersion and is typically of the order of 1000 km s 1 . From Eqs. (2) 
and (3), we find 

p fim p a 2 
Pe = Pg, P 



kT ( 



(5) 



- 



He nce the gas di stribution differs just by the index (3 in comparison with that of member galaxies. 

iKind (1 19621) derived an analytic approximation to the isothermal sphere of self-gravitational 
isothermal collision-less particles. The density profile of the cluster's member galaxies has been 
found to be well approximated by the King profile, 



Pg ~ PKing — Pa 



1+1- 

r r . 



-i -3/2 



(6) 



Here r c represents a core radius within which the density is regarded as constant. From Eqs. (4) and 
(5), the isothermal gas distribution is represented by: 



-| -3/3/2 



(7) 



This formula is called the isothermal /3-model (ICavaliere & Fusco-Femianol 19761) . 



2.2 X-ray emission process 

An X-ray spectrum emitted from an ionized intracluster plasma is described with a combination of 
continuum emission and line emission from heavy elements. The former is produced by free-free (or 
bremssstrahlung), free-bound, two-photon emission and the latter is by one -electron radiative tran- 
sitions, dielectric recombination satellite lines, and inner-shell ionization (Ivan Paradiis & Bleekerl 
ll999tlBohringer & Wernerll2010l) . 

In the temperature range of clusters (kT > 2 keV), the total emission is dominated by the 
free-free emission if the abundance of heavy elements does not exceed the solar abundance by very 
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Energy (keV) 

Fig. 2 APEC thin-thermal plasma model for fcT/[kcV] =1 (black), 2(red), 4(green), 
8(blue), 16(cyan), and the metal abundance of 0.3 solar. The graphs are shifted in the 
y-direction for clarity. 



much. The emissivity of the free-free emission at a frequency v from a hot plasma with an electron 
temperature of T g is given by 

^ = S$ ( 2^fc )1/2ne ^ Z2riifffr(Z ' r9 ' !y) X T 9 " 1/2ex P(-W^) (8) 

= A(T,Z,v)nl (9) 

where Z is an io n of charge in a plasma, and n e are the number density of ions and electrons, 
respectively (e.g. jRvbicki & Lightman|[l986l) . The Gaunt factor is a correction factor for quantum 
mechanical effects and is approximately gs ~ 0.9(hu/kT)~ 03 . The bolometric emissivity is then 

/>OG 

e ff = / e*dv = A(T, Z)n\ 
Jo 

~ 1.435 x 1CT 27 gTl ,2 n e ^ Z 2 n, [ergs" 1 cm" 3 ]. (10) 

Precise X-ray emis sion spectra from thin-thermal plas ma can be calculated by utilizing plasma 
codes such as APE C dSmith et al J 1200 lb . and MEKAL dMewe etalJll985l 1 19861: iKaastral 1 1 992t 



Lieda hl et al Jl995h . The updated version of the latter is available in the SPEX package dKaastra et alj 



19961) . For reference, the APEC models for various temperatures are plotted in Fig. [2] The metal 
abund ance is assumed to be .3 solar, as is typical of in tracluster gas dMushotzkv & LoewensteirJ 
Il997h . The abundance table of I Anders & Grevessd d 19891) is used her43. 

2.3 /3-proflle and hydrostatic mass estimate 

The surface brightness profile of an isothermal spherical plasma with a radial density profile given 
by Eq|7]is calculated by integrating the local emission per unit volume (EqlTOt and the density along 



The updated table for the solar system abundance is given in lLoddersI 120031) 
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the line of sight. We obtain the X-ray surface brightness S(r) at a projected radius r, 



S(r) = So 



S = n e0 n m A(T,Z) 



-3/3+1/2 



V¥r c T(3/3 - 1/2) 
AnDl T(3/3) 



[erg s 



(11) 



(12) 



Here n e o and uho are the central electron and hydrogen densities of the intracluster gas respectively 
and Dl is a luminosity distance to the object. It is known that the observed cluste r's X-ray surface 
brigh t ness is well fitted with the above function, and (3 ~ 0.6 — 0.7 on average (e.g. Jjones & Formanl 
11984 lOta & Mitsudall2004 . 

Once we have obtained the (3 profile parameters to characterize the surface brightness distribu- 
tion, we can estimate the three-dimensional density profile of the gas, i.e. the /3-model. Then the 
mass of the gas inside a radius r is given by integrating Eq|7] 



M &as {r) = 47rp gas (0)r c 



(1 



(13) 



where x = r/r c . From the assumption of hydrostatic equilibrium, we derive the total mass of the 
cluster inside a radius r, 

3kT/3r (r/r c ) 2 



M{r) = 



(14) 



pm p G 1 + (r/r c ) 2 

If gas is not isothermal and its temperature distribution has a radial dependence, T g (r), the hydro 



static mass is rewritten as 



M{r) = - 



kT g( r ) r 

G^JLTUp 



d In n g (r) d In T g (r) 



din r 



d In r 



(15) 



As an illustration, the result of hydrostatic mass estimation under the isothermal ,5-model for 
a gravitational lensing cluster CL0024+17 (z = 0.395) is shown in Fig. [3] and compared w ith an 
independent mass determination based on the gravitational lensing effect dTvson et al. 1 119981) . Sance 
the lensing effect directly maps the surface mass density of the cluster, regardless of the internal 
dynamical and thermal state of the cluster, comparison of the two methods provides information on 
the physical state of clusters dHattori et al.lll999t iKneib & Natarajanll201 ll for review). In the case 
of CL0024+17, a factor of two-three discrepancy has been found between the hydrostatic and strong 
lensing mass estimates, indicating that the system is experien cing a line-of-sight merger dOta et alj 



lensing 
2004; Z 



Zhang et al.ll200l Ijee et al.ll2007fc IZu Hone et alll2009l) 



2.4 Universal dark matter density profile 



Navarro et al. (1997) found from their numerical simulations of structure formations under the Cold 
Dark Matter (CDM) model, that the collapsed dark matter halos with masses over several orders of 
magnitude follow a universal density profile, 

Ps 



PNFw(r) 



(r/r s )(l + r/r s y 



(16) 



where p s and r a are the characteristic density and length, respectively. p s is related to the critical 
density of the universe p C rit and the characteristic density S c through p s = 5 c p C rit. Instead of the flat 
core of the King profile, the NFW profile has a core with oc r _1 dependence. Although the density 
diverges at the center, the mass inside a radius r, 



M NFW (r) = 4irp s r s 



ln(l 



1 + x 



x = r/r 



(17) 
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Fig. 3 HST/WFPC2 image of a lensing cluster with multiple lens arcs, CL0024+17 (left). 
Mass density pro files for the /3-model and the NFW-model of the CL0024+17 main cluster 
( Ota et al1l2004l) are shown and compared with the lensing mass model by iTvson et al.l 
il998j) (right). 



converges to as r — > 0. 

The density distribution of intrac luster gas in hydrosta tic equilibrium with the NFW dark matter 
potential was analytically derived bv lMakino et al. d 1998b assuming the masses of gas and galaxies 
are negligibly small compared to the dark matter. 



—B 1 - 



ln(l 



Pgas(^) = /OgasOexp 



Then the mass of the gas within a radius r is given by, 



_ AirGpm p p s r 2 s 
B = kf 



(18) 



M gas (r) = 47rp gas (0)e- 



-B r 3 



x 2 {l + x) B l x dx. 



(19) 



ISuto et all d 1998b generalized the universal density profile to the form p oc + x")] x 

and numerically computed the gas density profile in hydrostatic equilibrium for the case of p = a, 
v = 1, and A = 3 — a with the restriction 1 < a < 2. Note that the case with a = 1, v = 1, 
and A = 2 corresponds to the NFW model. They further computed the X-ray surface brightness 
distribution at a projected radius r on the sky, and derived a useful fitting formula in the following 
generalized shape. 



S{r) oc 

<i>c = 

c = 

v = 



i 



r s <\>c 
0.3(2/a- i; 
0.41 - 5.4(2 
-0.68 - 5.09(a - l) 2 + (0.202 + 0.0206a 8 )B 1 



a) 6 + (0.585 + 6A7a~ 51B )B- a6 / 30 



These are valid for 5 < B < 20 and 1.0 < a < 1.6 in the range 10 < 4> < 4>i 
S^max) = 10 _4 S'(0). We refer to the formula with a = 1 as the SSM model hereafter. 



(20) 



where 
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The SSM model has a surface brightness distribution similar to the j3 profile over a wide range 
of r, although it has an excess over the /3-profile i n the central region b ecause of the strong concen- 
tration of the dark matter halo of the NFW model. iMakino et aTl (Il998l) fitted simulated gas profiles 
which obey the universal dark matter profile with the /3-profile function, and noted that the best-fit 
relation between the scale parameter and the /3-model core radius is given by r c = 0.22r s . 

Two kinds of density profiles, the /?-model and the SSM model, have been introduced so far 
since they give reasonable approximations to observed gas profiles in studying the global cluster 
structure. Deviation from those models sometimes seen at the center of clusters will be mentioned 
later. 



2.5 Formation of clusters and the virial radius 

Numerical simulations based on the CDM model predict hierarchical structure formation, so rather 
conti nuous accretion of matter and sub-cluster merging occur in the process of cluster formation 
(e.g. jMoore et alj|200ll) . Hence clusters reside at junctions of cosmic filaments and are connected 
to the surrounding filamentary structures. It is, however, practically important to define a 'cluster' 
based on some simple model. In this section we b riefly review the collapse scenario according to 
the spherical collapse model dGunn & Gottlll972h . This model predicts a very important physical 
quantity of clusters, the virial radius. 

At some time epoch, a certain region of the Universe which happens to have a higher mass 
density than the background due to fluctuations starts breaking away from the general expansion, 
and eventually collapses to form a cluster of galaxies. Since at the collapse epoch, SI ~ 1, we can 
neglect the A term in the equation of motion of the shell. We also assume that the amplitude of the 
density perturbation is small, i.e S <C 1. Then we have 

d 2 r GM 

w = — (21) 



where M is the mass inside the shell and is constant. The first integral of this equation is, 

•\ 2 2GM 

- = +C. (22) 



dt I 



C is a constant, and the total energy C/2 must be negative for collapse to occur. The solution of 
Eq. d22b is given in a parametric form, 

GM ,„ . „ GM 

r = — r (l-cos0). (23) 



|qa/2v » \ C \ 

The radius, r, is at 8 = 0, i.e. t — 0. Then it increases with increasing 9 and takes the 
maximum,r m = 2GM/|C| at = w, i.e. ,t = t m = ttGM/|C| 3 / 2 (turn around). Then it shrinks 
to again at = 2ir, i.e., t = t c = 2irGM/\C\ 3 / 2 (collapse). After collapse, the system will be 
virialized. In the virialized system, the potential energy is related to the total energy as W = 2E. 
Assuming the radius of the system after virialization is r wlr , we have 

^ = -^ = 2^-2^!. (24) 

' vir 

Therefore, r v ; r = r m /2. The average density inside the virial radius r v - 1T is 

3 |Cf 



4tt G 3 Af 2 ' 



(25) 
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On the other hand, the solution of Eq. (1221 with C = describes the background expansion, 
because f2 ~ 1. The solution is 

q \ 1/3 , 

-GMJ t 2/3 . (26) 



The density inside rb gives the critical density at t, 



Pcrltit - Q - 24 ^ 3 qJ M 2- ( 2? ) 

Thus we obtain the important relation, 

A c = = 18tt 2 . (28) 

Pcrit 

We can assume that a cluster is virialized within the overdensity radius ta at which the average 
density is equal to A c times the critical density of the collapsed epoch. 

The spherical collap se in an il + A = 1 Universe is presented in the Appendix of 
iNak amura & Sutol dl997l) : a fitting formula for the overdensity in the flat Universe with finite A 

is 

A c ~ 187r 2 fi a437 . (29) 

By taking A c = 180 or 200, the overdensity radius of riso or T2oo is often quoted as a measure 
of the cluster's virial radius. r$oo is also frequently used for the reason that there i s an indication from 
numerical simulation that the hydrostatic assumption is valid within that radius dEvrard et all 19961) 
as well as that X-ray signals being detected out to rsoo or roughly ~ 1 Mpc in many clusters (beyond 
that deeper exposure is required to trace emission from the tenuous matter). The tempe rature scaling 
for th e overdensity radii for various A is derived using a nearby X-ray cluster sample dArnaud et alj 
120051) : for A c = 500, it resulted in r 500 h(z) = (1104 ± 13)(fcT/5 keV) a57±002 kpc. 

2.6 Radiative cooling of gas 

Since hot intracluster gas loses its thermal energy via X-ray emission, radiative cooling may affect 
the cluster structure once the gas is settled in the cluster's potential. 
The thermal energy loss is expressed by 

dE e ~ 

~W = " £ (30) 
where E e is the thermal energy of electrons per unit volume and E e = 3n e kT g /2. The volume 

emissivity, e ff can be denoted as e ff = qgn^Tg^ 2 . Thus if the hot gas cools, keeping the density 
constant, the temperature decreases according to the following equation 

dT 9 _ _ nT l/2 „ _ ^Sj} 

dt - 9 ' 3fc 

The solution is 



1/2 ^qsn e 
aTg , a ~ — — ■ ( 31 ) 



Tg (t)^=T g (Q)^-^t. (32) 



Thus the hot gas cools on the time scale 



fccool — — (JJ) 

a q s n e 

1/2 



sxiovf^V (- 



(34) 

The cooling time scale for the central region of typical relaxed clusters is estimated to be shorter 
than the Hubble time, the cluster core may be subject to radiative cooling. On the other hand, the 
radiative cooling is not considered to be important outside the core region because of lower density. 
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3 THERMAL EVOLUTION OF INTRACLUSTER GAS 
3.1 Cooling problem 

According to Section l2~6l gas at the cluster's center can radiate an amount of energy comparable to its 
total thermal energy in less than the Hubble time and thus cools. It was suggested from earlier works 
that the "cooling-flow" phenomenon would occur if the gas cools isobarically and no heating pro cess 
balances this cooling, so the gas flows inward maintaining the thermal pressure dFabianll 19941) . On 
the other hand, high-resolution X-ray spectroscopy showed that the temperature drops in the cooling 
cores by only a factor of two-three, and th ere is much less em ission at low temperature, as predicted 
by the standard cooling-flow model (e.g. jKaastra et al1l2004l) . This observational finding triggered 
explorations into a variety of scenarios for gas heating: heat conduction, active galactic nucleus 
(AGN) heating, magnetic reconnections, cosmic -ray heating etc. X-ray an d radio observations hav e 
provided evidence for the interaction of AGN jets with cluster gas (e.g., lMcNamaraetai1l2000l) . 
Although the work done by uplifting AGN bubbles on the surrounding gas may be of the order of 
magnitude to compensate the radiation loss, how the feedback achieves a tuning between cooling 
and heating is not clear. The similarity and smooth ness of cooling profiles in dicate the need for a 
continuous, distributed heat source (for review, e.g., |Peterson"&^ Fabian 2006). 



3.2 Statistical properties of cluster cores 

Regarding the density profile, a deviation from the conventional isothermal /3-model is commonly 
seen at the center of clusters having a compact core (often termed Cool Core (CC) clusters): 
they exhibit system atically higher central den sity while the profiles are fairly universal outside 
0-lr5oo ~ 100 kpc ( Neumann & Ar naud 1999). Fig.|4]shows the gas density profiles derived with 



the single-/? model (lOta et alj|2006 ). The density scatter is prominent within ~ Q.l 7- 500 and is found 



to be a significant source of scatter in the X-ray luminosity-temperature correlation ( Ota et al]l200d : 
lO'Hara et al.ll2006l: IChen et~ai1l2007l) . 

The statistical properties of gas density s tructure have been investigated from systematic anal- 
ysis of cluster samples by many authors (e g.JO'Hara et al .1120061 lOta et al]l2006t IChen et alJl2007t 



ISantos et al]|2008t ICavagnolo et alj|2009t iHudson et al.l l2010h . These X-ray studies show that the 



fraction of CC clusters is roughly 50% . The rest of the samp l e with out the central cool emission is 
called Non-Cool Core (NCC) clusters, lota & Mitsudal d2002[ 120041) first pointed out from the anal- 
ysis of ROSAT and ASCA archival data that the histogram of the cluster core radius exhibits a high 
concentration around 50 and 200 fey^k pc (Fig. [5]>. Later, a similar double-peaked distribution of 
core radius was shown independently bv lHudson et alj d2010l) : they utilized the Chandra data set on 
a nearby flux-limited sample with higher resolution. The consistency between the two results gives 
a confirmation of this nature. 

The relaxed clusters often host a central dominant elliptical galaxy, a lso called a cD g alaxy, 
which deepens the cluster potential well and causes a peaked gas profile (llkebe et alJ l 1999). The 
regular clusters with a small core tend to contain a cD galaxy, however, not all of them have one. 
Thus it is unlikely that the small core represents the potential distribution of the cD galaxy itself 
dAkahori & Masaill2005l) . 

Under the self-similar model, the internal structure of the gas should be scaled by the virial 
radius, and then r^/r c should be constant for all clusters. However, r c does not simply scale by 
r 500 (Fig. ID, particularly for those having a small core radius (r c < 100 kpc). This clear departure 
from the self-similar relation for small-core clusters suggests that the formation of the small cores is 
determined by some physical process other than the self-similar collapse. 

An investigation of X-ray fundamental plane gives another clue to explore the evolution of 
clusters. The presence of a planer distribution of nearby clusters in 3 -dimension al parameter space 
(the central gas density n r o, core radius r c , and temperature T) was first noted bv lFuiita & Takaharal 
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Fig. 4 (Left) Electron density profiles for 69 clusters. The best-fit density profile derived 
with the single /?-model are plotted, where the radius is normalized with rsoo- O.lrsoo 
is indicated with the vertical dotted line, inside which the scatter is the most prominent. 
(Right) Lx — T relation of clusters. A significant offset in the normalization factor of 
the Lx — T relatio n between clust ers with small (r c < 100 kpc) and large core radii 
(> 100 kpc) is seen dOta et alj|2006l) . 

Jl999b . implying t hat the clusters f orm a two-parameter family. Applying this technique to distant 
clusters at z > O.l. lOta et al.l 12006) obtained the following three orthogonal parameters: 

X oc n° 44 r°- 65 r-°- 62 (35) 
Y oc n°J 5 r° c A4 T - 78 (36) 
Z oc n° e J s r-°- 62 T- - 10 (37) 

and also confirmed the presence of the X-ray fundamental plane for the distant cluster sample. The 
distribution of clusters projected onto the X — Z plane is shown in Fig. [6] The Z-axis of the plane 
is called the principal axis and represents the direction along which the dispersion of the data points 
becomes the largest in the 3D space. By setting X ~ constant, Eqj37]yields Z oc r - 1 -78j'— o.io (x 
UgQ 20 ^ -0 ' 69 ■ Since the radiative cooling time is i coo i oc T x / 2 n~Q, it is rewritten as 

Z ex ■ (38) 

Therefore i coo i is considered to be a key parameter to control the cluster's gas evol ution. A trend o f 
morphological chan ge of X-ray clusters along the i coo i-axis is actually observed dOta et alj|2006l) . 
iHudson etal]d2010b noted that the cooling time is the most suitable parameter to segregate CC/NCC 
clusters, which is in agreement with the above result. 

It should be noted that a cool core is also found in some irregular clusters. This phenomenon is 
interpreted as a remnant of a merging core and may be used to diagnose the merging history (e.g., 
iMarkevitch & VikhUninll2007l) . 

3.3 Beyond the /3-modeI: observed gas profiles and the possibility of quasi-hydrostatic cooling 



To better reproduce observed X-ray surface br ightness profiles, some authors have introduced em- 
pirical models such as the double-/? model (I Jones & Formanlll984l) and the modified [3 model 
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Fig. 5 (a) r5oo — r c relation and (b) histogram of r c for 69 clusters at z > 0.1. In panel 
(a), 35 clusters with a small core of r c < 100 kpc and 34 clusters with a larger core 
of r c > 100 kpc are shown with the asterisks and the circles, respectively. The dotted 
lines indicate the self-similar condition corresponding to four different constant values of 

?"50o/ r c- 



(a) Distant sample 



* r<10()kpc 
O r>100kpc 




-2 
logX 



(b) Nearby sample 




-2 
logX 



Fig. 6 X-ray fundamental plane for the distan t cluster sample dOta et alj|2006l) and the 
nearby sample compiled bv lMohr et al.l dl999l) . The distribution of the clusters projected 
onto the l og X — log Z plane is shown in each panel. For the nearby sample, according to 
Table 2 o f lMohret al.l (1 19991) . non cooling-flow clusters are shown with the red triangles, 
and inner-core and outer-core components of cooling-flow clusters are separately shown 
with the solid blue boxes and open black boxes. 
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(|Vikhlinin et alJl200d) . For the latter, they modified the original /3-profile by adding a cool density 
cusp at the center and steepening the slope at a large radius. This model gave a good fit to spatially- 
resolved spectroscopic and imaging data t aken with Chandra for its large radial range including the 
core emission (see also lBulbuletal1l2010l) . 

Given that radiative cooling plays an important role in the thermodynamical evolution of ICM, 
how are the cool cores actually formed an d maintained? The possibil ity of quasi-hydrostatic cooling 
in the cluster core has been first noted 

by|M asai & Kitavamal (|2004|) . Unlike isobaric cooling flows 
that increase the local density so the thermal pressure P(r) counteracts the local cooling, quasi- 
hydrostatic cooling allows the gas to modify its profile or core size so VP(r) balances the gravi- 
tational force. The inflow is so moderate that the hydrostatic balance is not disturbed significantly. 
The quasi-hydrostatic model predicts a temperature profile that approaches a constant temperature of 
~ 1/3 that of ambi ent, non-cooling gas , which agrees with those derived from X-ray observations 
of relaxed cl usters p Caastra et al. 2004; lAllen et al.luOOlt iTamura et al.l 12.00 lb . Using a hydrody- 
namics code. lAkahori & Masail d2006l) investigated the evolution of the core structure of radiatively 
cooling gas. They suggested a radiative-cooling origin for the appearance of a small (r c ~ 50 kpc) 
core, while cooling is not important in clusters with large cores. Their simulations also showed that 
the cluster core maintains the quasi-hydrostatic condition before the initial central cooling time has 
elapsed. This result gives a possible interpretation on the observed double-peaked distribution of 
core size. 

lArnaud et all (1201 Oh discussed the universal pressure profile for the REXCES S cluster sample, 
and obtained the best-fit profile bas ed on the generalized NFW model bv lNagai et al.l (120071) . For the 
scaled temperature and the density. lArnaud et al.ld2010|) found that their deviations from the average 
scaled profile are anti-correlated with each other in the core, r/r^ < 0.2 (figure 3 in their paper); 
the anti-correlation is more clearly seen for cool core clusters. This behavior is supported by the 
quasi-hydrostatic cooling picture. 

Since the cooling time is shorter than the Hubble time for the CC clusters, some heating 
is needed to sustain the system, otherwise it would disappear ~Gyr after becoming virialized. 
Practically, however, heating due to mergers is likely invoked in the cluster's evolution. The clusters 
of core radii > 400 kpc in the histogram (Pig. [5j are attributed to mergers from their irregular mor- 
phology. Recently, th e process of cyclic evolution between CC and NCC clusters was proposed by 
iRossetti et al.l (|201 lb taking account of the lifetime of diffuse radio emission. 

3.4 Entropy profiles 

Measurement of a gas entropy profile provides important information on the evolution of gas since 
it determines the structure of intracluster gas and records the thermal history. The gas entropy, S, in 
the field of cluster research is defined by 

S = kTn- 2/3 , (39) 

and is different from the original definition in the field of thermodynamics. 

The gravitational heating, namely conversion of the potential energy to thermal energy, should 
depend on the depth of the gravitational potential, which is approximated by the virial temperature 
of the system. The entropy generation du e to gravitational collapse is pred icted to be self-similar 
and follow a power-law form S(r) oc r 1 1 dTozzi & Normanl200ll : IVoitll2005l) . Thus deviations from 
this baseline distribution may be attributed to cooling and heating processes in the cluster. Earlier 
results on groups and clusters observed with ROSAT showed that smaller systems like groups have 
entropy exc ess called the "entro py floor" at the center while the slope of the distribution follows the 
oc r 11 law JPonman et al.ll200"3l) . Thus the non-gravitational effects, preheating or galaxy feedbacks, 
are considered to play a greater role in smaller systems. 
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Fig. 7 Radial temperat ure and entropy profiles of the Hydra A cluster measured with 
Suzaku iSato et al.ll2012l) . Th e results for filam ent (black) and void (red) directions are 
shown. The Chandra results dDavid et al.l200lh are plotted in blue. 



Mor e systematic stud i es of entropy profiles with a large number of clusters have been car- 
ried out; ICavagnolo et al.l d2009l) derived radial entropy profiles of ICM for 239 clusters with the 
Chandra data (the ACCEPT sample) to find that most entropy profiles are well fitted by a model 
consisting of a power-law plus a constant, K(r) = K + ATi O o(?yi00 kpc) Q . The best-fit pa- 
rameters are (K , K 100 , a) = (16.1, 150, 1.20), (156, 107, 1.23) for clusters with K < 50 and 
A'o > 50 keVcm 2 , respectively. They also showed that the distribution of central entropy A' is 
bimodal, which peaks at A'o ~ 15 keV cm 2 and ~ 150 keV cm 2 . A si milar two peaked distribution 
has be en found in the REXCESS sample observed with XMM-Newton dPratt et al.ll2010h . |Pratt et alJ 
(120061) measured the entropy profile in relaxed clusters to find that outside 0.1? 2oo the scaled en- 
tropy profile is consistent with gravitational heating while the scatter increases with smaller radius 
and suggested that the results agree wit h models of accretion shock. 

The advent of the Suzaku satellite dMitsuda et al.ll2007l) enables the measurement of gas prop- 
erties out to large radii because of its low background level and high sensitivity. The temper- 
ature and entropy distributions up to the virial radius have been derived for massive clusters 



ature ana entropy distributions u p to tne virial radius nave b een deri ved tor massive clusters 
dGeoree et al.ll2009l iReiprich et all2009HBautz et al j2009UKaw aharada et al.ll2010b : Hoshino et al.1 



l2010HSimionescu et alj|201 lHAkamatsu et al.ll201 ll) . The latest result of Hydra A dSato et al.ll2012l) 
is shown in Fig. [7] For those clusters observed with Suzaku, a systematic drop in temperature by a 
factor of about three from outside the core to r2oo was found and the entropy profiles become flatter 
beyond r5oo in comparison with the r 11 profile. Some explanations for observed low entropy are 
proposed and discussed: in-fall ing matter retains some of its kinetic energy in the form of bulk mo- 
tion (?), a gas cl umping effect (Simionescu et al.ll2~01 ll) . and deviation of electron temperature from 
ion temperature dAkahori & Yoshikawall2010l) . 



4 CLUSTER MERGER AND SEARCH FOR NON-THERMAL PHENOMENA 

According to the standard scenario of cosmic structure formation, clusters are believed to have grown 
into their present shape via collisions and mergers of smaller groups and clusters. A cluster merger 
has a kinetic energy of the order of 



E 



1 



■(Mi + M 2 )v 2 



10 65 erg 



Mi + M 2 
10 15 M ra 



.3000 kms- 



(40) 
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where Mi and M% are masses of two objects, v is the collision speed and v = 3000 kms -1 corre- 
sponds to the mach number of 2-3 in the intracluster medium. This is the most energetic event in the 
Universe since the Big Bang. If two such objects collide with each other under their mutual gravita- 
tional attraction, a huge amount of energy may be released and a certain fraction is expected to heat 
the gas and generate non-thermal particles through shock waves, and induce bulk and turbulent gas 
motions. 

We can recognize signatures of merging in many ways. In X-rays, irregular morphology and the 
complex temperature structure of gas tell us that the system is disturbed due to the past mergers. The 
most p rominent shock feature has been detected in 1E0657-56 (the Bullet Cluster; iMarkevitch et al.l 
(120021) ). The bow shock propagates in front of a bullet-like gas and significant jumps in temperature 
and density have been found. The displacement between the peak positions of X-ray gas and dark 
matter distributions have been identified in merging systems such as the Bullet cluster and A2744 
dMerten et al.ll2.Ql H) . which provide an opportunity to constrain the self-interaction cross section of 
dark matter particles. 



4.1 Gas bulk motion and turbulence 

In the course of merging, different portions of the hot gas are predicted to collide with each other 
at a re lative speed of ~ a few x 1000 km s _1 , which will persist for several Gyrs after each merger 
event dNorman & Brvanll 19991) . If the gas has a large bulk velocity compared to its sound velocity, 
non-thermal pressure can no longer be neglected and has to be taken into account in estimating the 
cluster's mass. Suppose, for simplicity, that the gas is rigidly rotating with a circular velocity of 
ay(oc r), then the balance against the gravitational pull at radius r on the rotational equatorial plane 
then becomes, 

where 

and / is the fraction of gas that is rotating dOta et al.ll2007l) . Therefore the hydrostatic mass needs to 
be modified by a factor of (1 + f/3 r ) given the presence of kinetic gas motion. 

It is essential to constrain the gas motion through observations. The cluster gas contains a large 
amount of heavy elements such as iron, silicon, and oxygen etc. If the gas has a velocity along the 
line of sight, it produces Doppler shifts in emission lines from the heavy ions. The line shift due to 
the line-of-sight bulk velocity Vt UT b can be expressed as follows 

A£l „, k = £ « = , reV (_^_)(_^ =T) , 

where Eq denotes the rest-frame energy of the line emission. For example, a line shift due to the 
bulk velocity of 1000 kms -1 corresponds to a shift in the 6.7 keV Fe-K line energy by 22 eV. On 
the other hand, line broadenings due to turbulent and thermal motions are given by the following 
two equations: 

AE tmh = E a ^ = 6.7 oV ( r A (44) 



^6.7 keV / V300kms 



Because the thermal width is inversely related to the ion mass m, the contrast of turbulent broadening 
against the thermal one, A£^ tU rn/A£'th, becomes larger for larger m. Thus the Fe emission line is the 
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Fig. 8 6.7 keV Fe XXV line profiles for kT = 5 keV thermal gas convolved with the 
typical CCD and calorimeter detector responses, (i) thermal broadening only (black), (ii) 
thermal broadening + bulk motion Ubuik = 300 kins -1 (blue), (iii) thermal broadening + 
turbulent broadening v tU rb = 300 kms -1 (red). 



best-suited for velocity diagnostics in clusters. Fig. [8] shows the Fe emission line model convolved 
with instrumental responses assuming an X-ray CCD resolution of 130 eV and an X-ray micro- 
calorimeter resolution of 5 eV (FWHM). 

By measuring the energy shift with X-ray spectroscopy, one can directly probe the dynamical 
state of the gas. However, it is not easy since it requires not only a high energy resolution and a 
good sensitivity but also a precise instrumental energy-gain c alibration. Based on the careful assess- 
ment of positional gain variation of the Suzaku XIS detector (iKovama et al.ll2.Q07h , a tight constraint 
on the bulk veloci ty with an accur acy of 700 kms" 1 has been placed in the central region of the 
Centaurus cluster ( Ota et al1l2007l) . They placed the upper limit on the line-of-sight velocity dif- 
ference as 1400 kms" 1 . Hence, cluster mass estimation under a hydrostatic assumption is justified 
within a factor of about two-th ree. The Doppler- s hift measuremen t using t he Fe line has bee n carried 
out in several nearby clusters: ISato et al.l ( p008l): ISugawara et all (12009b; ISato et al.l d201 ll) derived 
the upper limit on the bulk velocity, and iDupke etaL ( 20071) : Dupke & Bregman (12006 ) reported 



possible detection of bulk gas flow. Recently, the significant bulk velocity of a sub cluster region 
relativ e to the main cluster, ~ 1500 kms -1 , has been detected in A2256 by Suzaku dTamura et al.l 
1201 ll) . 

The turbul ent motion has been pr obed by measuring a spatially-resolved gas pressure map in the 
Coma cluster (ISchuecker et al.ll2.Q04l) . The pressure fluctuation spectrum is found to be consistent 
with the Kolmogorov spectrum, yielding the lower limit of 10% of the total gas pressure in turbulent 
form. The turbulent line broadening has been constrained using the Reflection Grating Spectrometer 
(RGS ) on XMM-Newton in the central regions of ellipticals, groups and clusters (ISanders et al.l 
1201 lb . They placed a strong upper limit on the turbulent motion (< 200 kms -1 ) for several ob- 
jects while line broadening has been found in Klemola 44 and a weak signature in RX J1347-1 145. 

Theoretical expectation for line shifting and bro adening associated wi t h turbulence and bulk 
motions as well as their de t ectability a r e disc ussed by ISunvaev et al.l d2003l) : llnogamov & SunvaevI 
(l2003l) : lDolag et al.l d2005l ): lPawl et al.l d2005l) . 
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4.2 Non-thermal hard X-ray emission 

At radio wavelengths, s ynchrotron emissions e xtending over a Mpc scale have been discovered from 
more than 30 clusters dGiovannini et al ] ll999h . The existence of radio halo emission suggests that 
relativistic electrons are being accelerated in the intracluster space. Interestingly, there is a correla- 
tion between the radio synchrotron power (non-thermal, Pi. 4) and X-ray luminosity (thermal, Lx) 
for merging clusters while relaxed clusters w ithout a radio halo l ie in a region well separated from 
the merging clusters on the Pi. 4 — Lx plane dBrunetti et al. 2009]). It is s uggested that generat ion of 
high-energy particles is connected to the dynamical evolution of clusters dCassano et al.ll2O10h . 

In X-rays, the same population of high-energy electrons are thought to interact with 3K CMB 
photons and then generate non-thermal Inverse-Compton (IC) emission. The IC emission in excess 
of the thermal emission is then predicted to be seen in the hard X-ray band (>~ 10 keV) where the 
thermal emission normally diminishes because of the exponential cutoff in the continuum spectrum 
(Eq. |8). In addition, from the radio observation alone, we cannot separate the energy of magnetic 
fields from the energy of high-energy electrons. However, by comparing the radio and hard X-ray 
fluxes (S'gyn and Sic), the cluster's magnetic field is also estimated under the assumption that same 
population of relativistic electrons scatter off of CMB photons since the ratio S syn /Sic is equal to 
the ratio between the energy density of the magnetic field and the CMB 



S 8yn /Sic = Ub/Ucmb (46) 

dRvbicki & LightmanlfT986l) . U B = {B 2 /8ir) and U C mb = 4.2 x 10~ 13 (1 + z) 4 erg cm" 3 . 
The exact derivations of the synchrotron and IC emissions at a certain frequency are presented in 
iBlumenthal & Gouldldl970l) . 

Th e existence of non-therm al IC hard X-rays in the Com a cluster has been pointed out from 
RXTE dRephaeli & Gruberll2002l) and BeppoSAX observations dFusco-Femiano et alj|2004l) . Recent 



reports based on the broad band X-ray observations with Suzaku ( Wik et alj 12009 ) and Swift 
dWik et al.ll201 ll) did not find any significant non-thermal hard X-ray emission and the hard X- 
ray flux is reproduced by thermal models. This mismatch among several sat ellites is suggested to 
be re conciled if different sizes of field-of-views are taken into consideration dFusco-Femiano et alj 

l20iTb . 

Non-thermal hard X-ray emission has been constrained in about 10 bright clusters with Suzaku. 
The Hard X-ray Detector (HXD) on Suzaku has a field of view of 34' x 34' (FWHM) at ener- 
gies below 100 keV and has achieved the lowest background level dTakahashi et alj|2007l) . Fig. 
shows the hard X-ray spectrum of the hottest Abell cluster A2163 (z = 0.203) obtained with 
the Suzaku HXD. The additional power-law component does not significantly improve the fit and 
the observed hard X-ray spectrum is well explained by the multi-temperature thermal model, giv- 
ing the upper limit on the IC emissio n. This is consistent with the previous report by BeppoSAX 



Feretti et al.ld200ll). The HXD results dKitaguchi et al.ll2007t iFuiita et al.ll2008HKawano et alj|2009t 



Nakazawa et al. 120091 ISugawara et al.ll2009t IWik et all 120091: iKawaharada et al.ll2010al) are com- 
pared with those from other satellites, RXTE, BeppoSAX, and Swift in Fig. [10] There is no firm 
detection of the IC emission reported for these 10 objects from Suzaku. The cluster magnetic field 
obtained through the synchrotron-IC measurement (Eq.|46|) based on the Suzaku HXD observations 
is also plotted in the figure. Note that the estimation of magnetic field may be affected by the as- 
sumption of index p of the electron distribution, N(j) = Nqj~ p (7 is the Lorentz factor of the 
electron). 

The situation of non-thermal X-rays from clusters remains uncertain, and higher sensitivity in 
the high-energy range is required to further explore the physics of gas heating and particle accelera- 
tion in clusters. 
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Fig. 9 Suzaku HXD spectrum of the hot cluster of galaxies A2163. Significant emission 
is detected up to about ~ 50 keV. The spectral model consisting of multi-temperature 
thermal plasma (many thin dotted black lines) plus a non-thermal power-law component 
(blue) is indicated. 



4.3 Super-hot thermal gas in violent mergers 

To study thermal structure in clusters offers important perspectives in understanding the merging 
configuration and heating process of the cluster's gas. Merger shock and evolution of temperature 
structure for ions and electrons have been studied by numerical simulations. Given high sound ve- 
locity in intracluster medium, it does not seem easy for an in-falling sub cluster to acquire high Mach 
number to form strong shocks. On the other hand, there are some pieces of observational evidence 
for strongly heated gas that is likely to be generated by high-speed (> 2000 km s _1 ) collisions. 

The presence of extremely hot gas in the most X-ray lu minous cluster R X J1347-1 145 has been 
confirmed by the Suzaku broad band spectroscopy (Fig. [TT] Ota et al. 2008]). From the joint analysis 
of the Suzaku and Chandra data, the temperature of a hot clump (Fig.Q]right) is measured to be about 
25 keV, which is more than two times higher than the surrounding gas. This unexpectedly high- 
tempe rature gas has been pointed out previously by observations of the Sunyaev-Zel'dovich(SZ) 
effect dKomatsu et al.ll200lt iKitavama etal]|2004 . and the broad-band X-ray data have improved 
the accuracy by 3-fold. Importantly, the X-ray spectrum of this hot component is more accurately 
represented by a thermal emission model rather than a non-thermal power-law model. The results 
support a scenario that this cluster has experienced a rece nt violent merge r as the very hot gas is 
over-pressured and predicted to be short-lived (~ 0.5 Gyr) dTakizawall 19991) . It is also worth noting 
that the super-hot thermal gas significantly contributes to the hard X-ray flux, which needs to be 
precisely modeled in the search for non-thermal IC emission. Under the detailed multi-temperature 
modeling of thermal emission components, non-thermal IC emiss ion is not found to be sig nificant in 
the ha rd X-ray spectra obtained with Suzaku for RX J1347-1 145 dOta et alj|2008l) . Coma dWik et alj 
120091) . Abell 2163 (Ota et al. in prep.) and the Bullet cluster (Nagayoshi et al. in prep.). 

Is the very hot gas commonly seen in merging systems? In a nearby merging cluster 
A3667, a similar hot ( > 13 keV) thermal component is suggested from the Suzaku observations 
dNakazawa et al. 1 12009b . Including a shock-front cluster, the Bullet cluster, the Chandra and XMM- 
Newton te mperature maps show tha t some clusters contain very hot (> 10 keV) gas. The Chandra 
analysis by iMillion & Alien! d2009l) indicated that the hard excess can be attributed to non-thermal 
gas or quasi-thermal gas with kT > 20 keV. 
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Fig. 10 Non-thermal IC hard X-ray flux and cluster magnetic field in 10 clus t ers ob 
tained with S uzaku ( red). For S u zaku, th e results are quoted fr om Kitaguchi et al. (2007) 
Fuiita et all (120081): lOta etail d2008l): i Kawano et all (12009b: iNakazawa et alf i200% 



Sugawara et al. I d2009h: IWik et all d2009b: iKawaharada et al l hoiOal) Ota et al. in prep 
and Nagayoshi et al. in prep. For Swift (blue) jAiello et al.l d2009j. 120 1 0b : IWik et alJfeOl lb 
For RXTE (magenta) and BeppoSAX(b\ack), see lRephaeli et al.l (120081) : references therein 



The collision velocity necessary to explain such super-hot thermal gas due to strong shock heat- 
ing is high (~ 3000 — 4000 kms^ 1 ), which challenges the Lambda CDM model of cosmology 
dLee & Komatsull2010l) . 

5 FUTURE PROSPECTS 

X-ray spectroscopy and imaging observations bring us rich information on the nature of galaxy 
clusters, not only the baryonic content but also dark matter that governs the mass structure of the 
objects. Large-scale cluster surveys in various wavelengths are now on-going or planned, aiming to 
reveal the structural evolution in the Universe and obtaining more stringent limits on cosmological 
parameters. The baryonic mass fraction and cluster abundance as a function of redshift have been 
used to constrain the dark matter and dark energy densities as well as the dark energy equation of 
state. These measurements require precise mass estimates of large numbers of clusters, and thus 
understanding the physical state of intracluster gas to calibrate scatter and redshift evolution and 
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Fig. 11 High-resolu tion SZ effect map taken at the 45-m Nobeyama telescope (left; 
lKitavama~e t al. 200% and Suzaku broad- band spectra of the most X-ray luminous clus- 
ter RX J1347-1 145 (righ dOta et aD2008l) . The XIS data below 10 keV and the HXD data 
above 10 keV are shown with crosses. The step functions show the best-fit thermal model 
consisting of multi-temperature components for the ambient gas (many black lines) plus 
the very hot thermal gas (red line) identified in the South-East region of the cluster (see 
also Fig.[TJright). 



uncov er any bias in relationships between cluster mass and observables (e.g.. iMaiumdar & Mohrl 
120031) . 

Overall, clusters are regular objects, having positive correlations between global quantities (gas 
temperature, bolometric luminosity, gas mass etc) and the total mass derived either from X-ray ob- 
servations or a gravitational lensing effect. However, deviations from the self-similar expectations 
have been observed in terms of the power-law slopes and scatters around them. They are considered 
to have originated from non-gravitational effects like radiative cooling, feedback from galaxies, bulk 
and turbulent gas motions, magnetic field support etc. 

Among these issues, measurement of velocity structure to high accuracy is expected to be car- 
ried out by future high-resolution spectroscop y using an X-ray micr o-calorimeter. The ASTRO-H 
satellite is scheduled to be launched in 2014 dTakahashi et alj|2010l) and will play an critical role 
in revealing the dynamics of clusters. The Soft X-ray Spectrometer (SXS) onboard ASTRO-H is a 
non-dispersive spec trometer and enables high-resolution (5 eV) observations for both point sources 
and diffuse objects dMitsuda et alj|20l"oh . SXS will measure the kinetic gas motions to an accuracy 
of ~ 100 kms" 1 thr ough observations of line emissions. The Hard X-ray Imager on ASTRO-H 
dKokubun et alj|2010l) will constrain the non-thermal hi gh-energy contents in clusters with its imag- 
ing spectroscopy in the hard X-ray band. Now NuStar dHarrison et alj|201oh is successfully in orbit 
and draws peoples' attention to upcoming observations with the first focusing telescope in the high 
energy X-ray regime. The eROSITA on the Spectrum-Roentgen-Gam ma mission will pe rform an 
all-sky survey in the X-ray energy range and detect ~ 100000 clusters dPredehl et al.ll2O10h . In con- 
junction with optical and SZ surveys, the next-generation X-ray missions will largely enhance the 
study of clusters and lead us to draw a more complete view of structure formation and evolution in 
the Universe. 
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